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The most compelling and popular models for dark matter predict that it should congregate and 
annihilate in stellar cores. Stars where annihilation contributes substantially to the total energy 
budget look very different to those with which we are familiar Here I explain the general features 
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^ of stars modified by dark matter annihilation with the help of a series of grids of 'dark' stellar 

^ evolutionary models, and describe the public code with which they were computed. I go on 

'— ' to discuss possible impacts of dark stars on the high-redshift Universe, including the history of 

reionisation. The preUminary reionisation calculations reproduced here are based on dedicated 
models for dark star atmospheres, and for the stellar populations to which dark stars would belong. 
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1. Introduction 

Weakly-interacting massive particles (WEVIPs) are amongst the most promising candidates for 
dark matter [1, 2, 3]. In the standard scenario, WIMPs are Majorana particles that self-annihilate 
with a cross-section of ~ 3 x 10^^^ cm^ s^^ leading to a thermal relic dark matter abundance very 
similar to the observed value. If such WIMPs are somehow confined within stellar cores, the 
resultant annihilation can have a substantial effect on stellar properties [4, 5]. The idea that such 
'dark stars' could exist has been pursued rather vigorously in recent years, with attention given to 
white dwarfs [6, 7, 8, 9], neutron stars [7, 10, 11] and main-sequence stars [12, 13, 14, 15, 16] 
in our own Galaxy, as well as to the first stars (Pop III) [17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 
27, 28, 29, 30, 31, 32]. The effect of dark matter annihilation on Pop III stars has been examined 
in the context of broader astronomical observations, including reionisation [33, 34], HST/JWST 
deep-field magnitudes and galaxy spectra [35, 36, 37]. 

For the purposes of these proceedings, and more generally, I define a dark star as 'any stellar 
object whose structure or evolution has been affected by dark matter annihilation'. Other authors 
prefer to use the term 'dark star' to refer exclusively to the dark-matter dominated phase that may 
occur during the formation of Pop III stars (e.g. [17]). The distinction is essentially just semantic, 
but important for readers' understanding of the literature, should they choose to delve into it. These 
proceedings are not intended as a comprehensive review of dark stars as a field of active research, 
rather as simply a set of summarising notes of the (semi-review) presentation given at CRF2010. 

Dark matter can find its way into stellar cores via two distinct physical mechanisms, as out- 
Uned in Fig. 1. The first is by the gravitational contraction of a baryonic gas cloud, as it cools 
and collapses during star formation. The baryonic contraction steepens the gravitational potential, 
drawing dark matter into the centre of collapsing cloud. The contraction of the dark matter need 
not be entirely adiabatic for this to occur, although adiabaticity is typically assumed. This mecha- 
nism can only occur during formation of a star, so the resulting population of dark matter cannot 
be replaced by the same mechanism after it annihilates away. The second mechanism is by nucleai^ 
scattering, whereby WIMPs scatter weakly off nucleons in the star, lose kinetic energy and become 
gravitationally bound to the star. They then follow bound orbits to return and scatter repeatedly, 
losing more energy and settling down to the stellar core. This mechanism can continue to be ef- 
ficient for as long as an appropriate population of dai^k matter exists for the star to capture from 
in the halo surrounding it; in this way, capture via nuclear scattering may in principle replenish a 
star's population of dark matter indefinitely. 

2. Structure and evolution of dark stars 

Regardless of the physical mechanism bringing dark matter into a star, the total amount 
present, and therefore the rate at which it injects energy into the centre of the star, are what ul- 
timately determine the impact upon stellai^ structure and evolution. Around 90% of the rest mass 
of each WIMP is injected into the gas locally as additional heat [14], so WIMP annihilation acts 
as an additional luminosity alongside nuclear burning. This additional term is implemented in 
the open source dark stellar evolution code DarkStars [16], along with highly detailed treatments 
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Figure 1: A schematic representation of the two physical mechanisms by which dark stars acquire dark 
matter. Gravitational contraction (left) occurs where baryonic collapse of a star-forming gas cloud steepens 
the gravitational potential, drawing dark matter in along with the baryons. Nuclear scattering (right) relies 
on weak scattering events between WIMPs and nucleons, whereby dark matter particles lose energy and 
become gravitationally bound to the star. 

of the capture of WIMPs, their distribution within the star, and conductive heat transport due to 
additional scattering events between bound WIMPs and nucleons in the the stellar core. 

Owing to the the negative specific heat of a self-gravitating system, the additional energy 
means that the cores of dark stars are cooler and less dense than those of normal main sequence 
stars. Dark stars are hence larger, cooler and more diffuse than main-sequence stars of similar 
masses, and in fact resemble protostars lying somewhere on the protostellar cooling (Hayashi) 
track. This resemblance can be understood in terms of the differing microscopic properties of 
the two energy sources. To a good approximation, the luminosity from WIMP annihilation is 
decoupled from the gas pressure or temperature, as the capture and annihilation rates are only 
weakly dependent upon the stellar structure. The same is true of gravitational potential energy, the 
primary energy source during protostellar cooling. Nuclear burning, on the other hand, is strongly 
dependent upon the core temperature and density. As WIMP annihilation becomes stronger, nuclear 
burning becomes progressively less efficient; the more dark matter is present in a stellar core, the 
more it expands and cools, and the further 'back up' the Hayashi track a star can be found in 
the HR diagram. This causes an overall decrease in the stellar luminosity for moderate WIMP 
capture rates, but an increase at large capture rates, where the annihilation energy alone exceeds the 
equivalent main-sequence luminosity. The partial replacement of fusion power with energy from 
dark matter annihilation also slows the rate at which hydrogen is converted to helium, increasing 
stars' main sequence lifetimes. This picture is summarised in Fig. 2, where we give evolutionary 
tracks for some example main-sequence stars following their immersion in dark matter halos of 
varying densities. 

Whilst the direct physical impacts of dark matter annihilation in stellai^ cores show only a very 
weak dependence upon metallicity [14], the fact that nucleai^ scattering and gravitational contrac- 
tion are effective at such different times in a stai^'s life cycle means that their impacts do end up 
being rather functionally distinct. In particular, the presence of large amounts of dark matter during 
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Figure 2: Example evolutionary tracks (left) and central equation-of-state diagrams (right) for a grid of 
low-mass Pop II stars, when evolved from the main sequence with different amounts of dark matter capture 
by nuclear scattering. From top to bottom, models exhibit sufficient capture to provide 10^^, 10^', 1 and 10 
times as much energy from WIMP annihilation as from nuclear fusion. For details, see Ref. [14]. Models 
were computed using the DarkStars code [16]. Qualitatively similar changes occur for higher masses and 
other metallicities. Stars that are strongly affected by dark matter annihilation before reaching the main 
sequence approach similar Hayashi-like solutions, but descend from higher on the Hayashi track instead of 
ascending to them from the main sequence. 
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the formation of a star (i.e. as can only be produced due to gravitational contraction) essentially 
inhibits the collapse, as it can partially prevent further cooling of the gas beyond a certain critical 
point, well before a main-sequence object has formed [17, 25]. This means that stars in which grav- 
itational contraction is an effective means of obtaining WIMPs take longer to halt gas accretion by 
radiative feedback, leading them to acquire larger masses than stai^s formed without the influence 
of dark matter. It is thought that because Pop III stars form near the centre of (almost) pristine 
dark matter halos, this effect should be essentially unique to Pop III stars. More speculative work 
[36] has even suggested that such a mechanism could lead to the formation of multi-million solar 
mass 'supermassive dark stars'. In any case, if Pop III dark stars do exist, their larger masses, in 
combination with the fact that the annihilation luminosity can exceed the equivalent main-sequence 
luminosity anyway, mean that they are not only cooler and larger than standard Pop III stars of the 
same mass, but also substantially more luminous. 



3. Impacts at high redshift 

If dark stars exist in the early Universe, they may be observable via their indirect impacts upon 
cosmological processes, or even directly as extremely red, highly luminous objects. As single 
objects, standard Pop III dark stars would only be observable with e,g. JWST if they were very 
cool and long-lived, and even then only when viewed through the strongest known gravitational 
lens [35]. Their unique spectral impact on integrated galaxy spectra could however provide a 
more promising detection channel [37]. Supermassive dark stars, on the other hand, should be 
individually observable by both HST and JWST in many scenarios [37], severely constraining their 
abundance and lifetimes. 

Because they can be more massive (and hence more luminous) than normal Pop III stars, and 
potentially live for far longer. Pop III dark stars could substantially change the reionisation history 
of the Universe. Exactly how reionisation is changed depends not only on the abundance and 
lifetimes of the dark stars, but to what extent they actually remain affected by WIMP annihilation 
over the course of the lifetimes. This is because large scattering-mediated WIMP capture rates 
lead to very cool dark stars, which do not produce a substantial number of ionising photons for 
most of their lives. This causes reionisation to be delayed relative to standard Pop III scenarios, 
as shown in Fig. 3 [34]. Here the evolution of the ionised hydrogen gas fraction in the Universe 
is shown as a function of redshift for a number of different stellar population models, with and 
without very cool dark stars. The presence of such cool dark stars can be seen to delay reionisation 
by different amounts, depending upon the dark star lifetimes and whether they are followed by a 
'normal' (non-dark) Pop III starburst, or simply by standard Pop II star formation. 

On the other hand, in the absence of large scattering-mediated capture rates, the larger stellar 
masses of dark stars compared to normal Pop III mean that they produce more ionising radiation 
than would otherwise be expected. This would hence have the opposite effect on reionisation, 
causing the Universe to become fully ionised at much higher redshifts than in the canonical Pop III 
scenario [33]. Studying the integrated optical depth of the CMB with WMAP and Planck should 
thus help to further constrain the abundance and lifetimes of the Pop III dark stars [34]. 
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Figure 3: Reionisation histories (in terms of the Hll fraction) from Ref. [34], computed under different 
assumptions about the populations of dark stars existent in the early Universe. Dark star curves include an 
initial population of highly dark-matter-dominated Pop 111 dark stars, which do not contribute any reionising 
photons, and hence delay reionisation. This population is assumed to live for either f = 10 or 100 Myr 
Immediately before their deaths, the dark stars briefly resemble ~8OOM0 main-sequence Pop 111 stars; 
atmospheres and reionising photon fluxes of dark stars at this stage are modelled using the TLUSTY code 
[38]. After the population of dark stars, reionisation models assume either standard Pop 11 star formation, or 
normal Pop 111 star formation followed by normal Pop 11 star formation. The normal Pop 111 star formation 
is modelled as having a duration of lOMyr, and consisting of stars with masses of 1OO-3OOM0. The 'no 
DS' model consists of just this normal Pop 111 for a duration of 10 Myr, followed by a normal Pop 11. The 
collapsed baryon fraction in halos is also shown for comparison. 
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